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We present a study of an Hanbury Brown and Twiss (HBT) interferometer realized with anyons. 
Such a device can directly probe entanglement and fractional statistics of initially uncorrelated 
particles. We calculate HBT cross-correlations of Abelian Laughlin anyons. The correlations we 
calculate exhibit partial bunching similar to bosons, indicating a substantial statistical transmuta- 
tion from the underlying electronic degrees of freedom. We also find qualitative differences between 
the anyonic signal and the corresponding bosonic or fermionic signals, indicating that anyons cannot 
be simply thought as intermediate between bosons and fermions. 



Two-particle interference is a major pillar of quan- 
tum mechanics, very much like the phenomenon of sin- 
gle particle interference. Such interference has been ob- 
served with photons in the historical Hanbury Brown 
and Twiss (HBT) experiment [1112], and much later with 
electrons [3]. Quantum Hall systems can exhibit emer- 
gent particles (dubbed anyons) with fractional statis- 
tics [H ini ■ Despite intensive study, direct signatures of 
anyonic statistics remain elusive. Here we study an HBT 
interferometer with anyons, which can directly probe en- 
tanglement and fractional statistics of initially uncor- 
related particles. Specifically, we calculate HBT cross- 
correlations of Abelian Laughlin anyons. The correla- 
tions exhibit partial bunching similar to bosons, indicat- 
ing a substantial statistical transmutation from the un- 
derlying electronic degrees of freedom [B]. Furthermore, 
we find qualitative differences between the anyonic sig- 
nal and the corresponding bosonic or fermionic signals, 
indicating that anyons cannot be simply thought as in- 
termediate between bosons and fermions. 

Edge channels of a fractional quantum Hall system of- 
fer a natural framework to study transport properties 
of anyons. Earlier attempts to consider entanglement of 
such quasiparticles (QPs) either addressed time-resolved 
correlation functions [7] (which may be very hard to mea- 
sure) or relied on a single source geometry setup [BHTO] 
(which may introduce superfluous interaction-induced 
correlations). Here we study zero frequency current- 
current correlations in a truly HBT interferometer setup, 
whose physics is governed by QPs dynamics. Because of 
their fractional charge and fractional statistics, scattering 
of these QPs results in non trivial correlations. Below, 
we consider the case v ^ l/i for concreteness, but, our 
analysis can be generalized to other Laughlin fractions. 

Consider first a heuristic estimate of these correlations, 
outlined in Fig. [T] Two particles are emitted respec- 
tively from two sources Si and 5*2 and scattered towards 
two detectors Di and D2 by a beam splitter, e.g. a 



quantum point contact (QPC) for electrons and QPs, 
or a half silvered mirror for photons. We evaluate the 
probability P(m, 2 — to), m = 0,1,2, that to particles 
are collected at the drain Di while (2 — to) are col- 
lected at the drain D2. Consider, e.g., the diagrams 
contributing to P( 1,1) [cf. Fig. [ija)]. Each diagram 
represents an amplitude contributing to P(l,l). Their 
weights are tt'exp[i(l/2)7ri/] and rr'eiq^[i{'i /2)ttv], with 
V — 0,1,1/3 for bosons, fermions, and v — l/i anyons, 
respectively. Note that we have included quantum statis- 
tics factors which reflect the extent by which one parti- 
cle winds around the other. It follows that P(l, 1) — 
|tt'exp[i(l/2)7rj/] -f rr'exp[i(3/2)7rzy] I . For simplicity we 
consider symmetric scatterers, |rp = |tp = 1/2, in which 
case P(l,l) = (1/2)(1 - cosvri/). Similarly, P(2,0) = 
P(0, 2) — (1/4)(1 -I- cosTTi/). For classical particles one 
sums up probabilities, rather than amplitudes, leading to 
Pc(l,l) = 1/2, Pc(2,0) = Pc(0,2) - 1/4. The results 
for fermions and bosons coincide with calculations based 
on second quantization |llj . For ^' = 1/3 this results in 
boson-like bunching [7] (P^=i/3(2,0) > Pc7(2,0)). 

Our main analysis, outlined below, reinforces the ob- 
servation that the scattering of two Laughlin anyons is 
bosonic-like. At the same time, it also reveals the non- 
analytic structure of the interferometry signal of such 
anyons, implying that the latter are not simple interpola- 
tion between fermions and bosons. The schematic setup 
is depicted in Fig. |2] What replaces optical beams in the 
solid state device are edge states of the quantum Hall ef- 
fect, formed due to the presence of strong perpendicular 
magnetic field. The chirality of these edge states allows 
the transport of charge excitations over large distances. 
The presence of an Aharonov-Bohm (AB) flux (^ab) pro- 
vides an important handle to control and analyze the 
HBT correlations. The analogue of half-silvered mirrors 
are quantum point contacts (QPCs), which facilitate con- 
trolled forward transmissioir/backscattering reflectioir of 
the impiirging charge excitatioir. Electroir interferome- 
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FIG. 1. Two-particle amplitudes contributing to: (a) P(l,l), 
the two particles are emitted from 5*1 and S2, and collected 
at Di and D2; (b) P(0,2), both particles are collected at Di. 
t,r,t',r' are the single particle scattering amplitudes, \t\'^ + 
|rp = l,|t'|^-|-|r'p = l. Note the statistical factors, reflected 
by the winding of one particle around the other. For (a) they 
are exp[i(l/2)7ri^] and exp[i(3/2)7r//], respectively. 

ters with such features have been realized [31 [T^]. We 
focus on the magnetic flux sensitive part of the current- 
current correlation, and show that the result is radically 
different from what has been predicted [13], and later 
observed [5] for the electronic case {1^ = 1). 

A QP in a quantum Hall liquid at Laughlin filling fac- 
tor v can be described as a composite object, consisting 
of a point charge q = ve with a single quantum mag- 
netic flux solenoid, $0 = hc/e, attached to it. When 
a QP encircles another QP it will pick up an AB-phase 
9 = 27rz/ which accounts for their mutual fractional statis- 
tics [13]. When a QP tunnels from the external to the 
internal edges, its flux is trapped inside the interferom- 
eter [151 US]- The magnetic flux enclosed in the active 
area of the interferometer (depicted in blue in Fig. [2| is 
*tot("') = ^•ab + ^statin), where <^stat is the statistical 
flux and is given by $0 times the number, n, of trapped 
QPs. The dynamics of QPs moving along the edges of the 
interferometer is then entirely determined by n mod(3). 
i.e. for a given value of ^ab the system can be found in 
three possible states characterized by n — 0, 1, 2. 

For the study of the non- equilibrium dynamics of our 
strongly interacting HBT interferometer, we address the 
Markovian evolution of the system among the three possi- 
ble values of the statistical flux. Our microscopic Keldysh 
analysis simplifies, and can be cast in terms of rate equa- 
tions for a certain parameter range [T7]. The rate equa- 
tions (whose coefficients are obtained by a microscopic 
analysis) carry information on interference effects of cur- 
rent cross correlations. Below we treat the QP tunneling 




FIG. 2. Scheme of a Hanbury Brown and Twiss setup real- 
ized with edge states. The relevant edges for our analysis are 
the lines SiDi with i=l,2,3,4. The external edges, S\D\ and 
S4D4, are kept at potential V, the internal ones, S2D2 and 
S3D3, are grounded. Inter-edge tunneling (dashed lines) takes 
place at the four QPCs, .4, B, C, T>, with tunneling ampli- 
tudes, r^, Fg, Fc, Fe), respectively. The bias V is assumed 
to be large compared to the thermal energy, fcsT, so that 
only processes that involve transfer of QPs from the external 
to the internal edges are relevant for our analysis. Li (not 
shown in the figure) is the distance between two consecutive 
QPCs along edge SiDi. The magnetic flux threading the blue 
area, "I'tot, is relevant for the interferometry discussed here. 



current at each QPC perturbatively. 

Let us define the quantities needed in the ensuing anal- 
ysis: {li) is the average tunneling current measured in 
drain i and S,.,, = J^, dt{{h{t) - (/,)) (/,,(0) - (/,.))) is 
the zero-frequency current-current correlations between 
drains i and i' . The latter is the main object of this Let- 
ter. Next, we define P{f,t\j), the probability to find the 
system with statistical flux / [indices are defined mod(3)] 
at time t given that it had statistical flux j at time zero. 
The system's dynamics is governed by a standard Master 
equation 

E [Pik,t\j)Wkj - PiLt\j)Wf.,k] . 

fc=0,l,2 

(1) 

Here Wjj is the total transition rate from the state 
J to the state /. In order to study the magnetic flux 
dependent part of the current-current correlations, we 
need to consistently include at least single-QP processes 
and two-QP processes, i.e. second and fourth order in the 
tunneling amplitudes F, respectively. In the limit of high 
voltage bias, eV ksT, considered here, only processes 
that transfer QPs from the outer to the inner edges are 
relevant. 

Several microscopic processes, labeled by C, contribute 
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to each Wjj such that Wjj = J2( ^jf- "^^^ processes 

allowed are either Wjj, wj^^+i, or wj^^^j- The former 
renormalizes the vacuum current and does not affect any 
quantity calculated below. j+i has contributions from 
single QP processes (independent of flux, hence, indepen- 
dent of j), as well as from two-QPs processes (dependent 
of flux). consists of two-QPs processes, and may 

or may not be flux-dependent. Each of the rates can be 
written as W^^j = W^'^^'k^^^^ with n^p discussed in the 
caption of Table |l] which depicts all relevant processes. 
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TABLE I. Elementary QP transfer processes. Each process, 
{Q — {m,N,(j)), is characterized according to the change, 
m, in the number of QPs trapped in the interferometer; 

the QPCs at which QP tunneling takes place; and the 
flux, 0, entering the flux factor k^'"'^''''-' = cos[27r(^/(3"l>o)]- 
Note that cf) = Q depicts a flux-independent process and 
4) = $tot(i) = 'I'AB + j ■ ^0 a process that depends on the 
total trapped flux. The order of the process (second or fourth 
in the tunneling amplitude F), the initial and final fluxon 
states [ij, /), where f — j is the added number of statistical 
fluxons], and the charge added at each drain {+1 refers to 
the absorption of one QP or charge q = — (l/3)|e| at the 
drain), are indicated. For example (cf. Fig. [3|, the pro- 
cess (C) = (1,^,0) corresponds to the emission of a QP 
from source Si, its tunneling across QPC A, and its trap- 
ping at Dz- Following the tunneling event a quasi-hole is 
created at edge S\D\ and a charge —q is consequently ab- 
sorbed in D\. The flux dependent processes [the two-QPs 
trapping process (2, ^Z3CI?, $tot(j)) and the single-QP trap- 
ping (l,^Z3CD,$tot(j))i] are illustrated in Fig.jsj 

Consider, first, the current collected at any of the 
drains. Assuming short-range interactions |18j . which is 
reasonable in the presence of a metallic top gate, this cur- 
rent is flux independent (similarly to the v = \ case |13)). 
hence, is not of interest for us here. The following argu- 
ment can be used to show this: consider for instance 




FIG. 3. (a) In process (2, y^lSCO, $tot(j)) two QPs are trans- 
ferred from edges 1 and 4 to edges 2 and 3, the process 
is AB-sensitive due to the interference between two ampli- 
tudes A\ and A-^,. In A\ a QP tunnels from edge 1 to 
edge 3 and a second QP tunnels from edge 4 to edge 2 
(red dotted line). In A-i a QP tunnels from edge 1 to edge 
2 and a second QP tunnels from edge 4 to edge 3 (blue 
dashed line). This process changes the statistical flux by 
two. (b) Process (1, ^SCO, $tot(i))i (and similarly process 
(1, ABCT), "l?tot(j))2) is also AB-sensitive but in this case only 
one QP is trapped inside the interferometer changing the sta- 
tistical flux by one. 

the current at drain D3: owing to the chiral propagation 
along the edges this tunneling current does not depend 
on the scattering at QPC T). A gauge transformation can 
then ascribe the total magnetic flux to QPC T) — hence 
the current in drain 3 is AB independent. A similar ar- 
gument holds for the tunneling currents collected at the 
other drains. 

We, next, consider the AB-dependent component of 
the cross-current correlations. It is sufficient to ex- 
press the following rates: the sing le-QP rates wf^^^Y = 
7|r_4p (and similar expressions for rates involving the 
processes (1,;B, 0), (1,C,0), and (1,1?, 0)); and the two- 
QPs rate 

^(2,^BCP,*.„.o)) ^ f7|r^r«rcrp|cos[^($AB+j-<i>o)/$o]. 

Here the P's are the QPs tunneling amplitudes at the four 
QPCs, 7 and fl are coefficients to be calculated below. 
Using the method developed in references [ini HO] we are 
able to calculate the AB-dependent component of the 
cross-current correlator ^i^: 

^1^ = ^lA ~ ('S'i,4)ab (2) 
^ e^\^A^B^c^v?^'' cos[27r($AB/<&o)] 

6(|r^|2 + |rB|2-H|rcP-i-|rp|2)272 ' 

where ()ab refers to averaging over <i>AB- 

Model and Methods — The low energy physics of the 
system is well described by the effective bosonic Hamil- 
tonian [l8] 

n * r 

Ho = ^Y.h<^-'^i)^^ (3) 

7—1 
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describing chiral plasmonic excitations on the four 
edges (e.g. SiDi) of the interferometer. Here, 
V is the plasmonic velocity at the edge. The 
bosonic fields satisfy the commutation relations are 
[4>iix,t = 0),(j)kix',t = 0)] = iTTSikSgn{x - x'). The oper- 
ators exp{i(j)i / y/u) and exp(i-\/i/0/) are, respectively, pro- 
portional to the electron and the quasiparticle creation 
operator on the edge I. 

To fully account for the quantum statistics of such par- 
ticles one needs to multiply these bosonic operators by 
"string operators" , known as Klein factors [3 [13 IHl HI] ■ 
In our analysis this procedure is replaced by carefully ac- 
counting for the dynamics of the statistical flux, attached 
to the tunneling QPs. 

The total Hamiltonian, H = Hq + Hx, includes a 
tunneling part, Ht = (H^ + + + H^) + h.c, 
which accounts for the most relevant tunneling opera- 
tors at the QPCs. We assume that the external (inter- 
nal) edges are tuned at voltage V (0), select a gauge 
whereby the flux dependence is attached to -ff^, and re- 
define the vacuum value of the fields at the external edges 
4>i{x) —i' (j>i{x) - eV ^x / {vh) (/ = 1,4). With these ma- 
nipulations the tunneling operators read 



Note that the magnetic flux attached to comprises of 
both the AB-flux and the statistical flux due to n mod{3) 
QPs. 

We next calculate the transition rates. The above 
model facilitates the calculation of the rates of the pro- 
cesses appearing in Table llj Rates are computed us- 
ing generalized Fermi's golden rule (see, e.g. [22 ) in 
order to evaluate single and two particles transfer be- 
tween the edges. Generally we can write the transition 
rate between any initial state {ipi) with thermal occu- 
pation Pi to any final state {ipf) obtained from the ini- 
tial one by transferring one or two QPs as W^/fj?/ = 
{271 /h)j:^^p,\{i;,\f\i;f)\^6{Ef - E,), where f ^ Ht + 
Ht{E.i ~Ho-iO+)-'^HT-\ . For example, 



f 

J — c 



r-,r{2,ABCT},f>tat (j ) ) 

dt 



2tt 

= |r^rBrcri,|cos[-^(<E'AB + <E>j)/$o]x 



"^Gt (e, -U)G'^ {€+ueV, U)G^ {e+ueV, L2)Gf (e, -Li) 



tion. In time-space representation they are given by 

G>{t,x) = (e-iv^0.(^.t)e»v^0,(o,o)^ ^^^^ Gf{t,x) = 



(4) 



Here we have introduced G<^{e,x) and Gf(e,x), 
the Green's functions in energy-space representa- 



^g^v/;70,(o,o)g-^^/?0,(^,t)). We find for 7 and 1) (cf. Eq. (|2 



7 
and 



C 



^l/3g7rQ/2 



[l + 2cosh(7ra)]r(f)r(; 



Or(| + f)' 



327ri/3^^/'/35/3 

2l/3r(2/3)a4/3;i7/3„4/3 



(5) 



(6) 



where a = veV/{-KkBT), /? = l/(fcsr), and C = 
{lc/v)'^^^(2T:/h)^/^. In order to obtain Eq.([6]) we have 
assumed that the lengths of the individual interferom- 
eter arms satisfy Li + L4 — L2 — L3 <^ hv/{i>eV) and 
a > 1. 

Discussion — Eq. ^ is our main result. We first note 
that the leading AB cross-current dependence comes with 
a plus sign, akin to bosonic HBT correlations (cf. refs. [7] 
and|23|. This conclusion is in qualitative agreement with 
our toy model discussed above. The structure of Eq. ^ is 
worth noting as well. For an electronic two-particle inter- 
ferometer operating in the integer quantum Hall regime 
(j/ = 1), the leading flux-dependent contribution in the 
weak tunneling regime is proportional to [13] F** (bosons 
would behave the same way). Likewise one might expect 
the fundamental flux periodicity of S^f to be $0, in line 
with gauge invariance [HI dni [H] . This would suggest 
that is proportional to F^^, representing a coherent 
sequence of three {2, ABCVj^totij)) two-QP processes. 

implying that an ex- 



Our result for scales as F 

pansion of S^^ in two-QPs rates, W. 



{2,ABCV,<S>tatU)) 



non-analytic. This unique scaling with F is the signature 
of QP HBT interference. Formally, this intriguing behav- 
ior is the outcome of the dressing of two-QPs processes 
by an infinite series of single-QP processes. We notice 
that the above results applies to the case of = 1/3 con- 
sidered here; for a generic Laughlin filling factor v = I /m 
(to odd), having m possible values of the statistical flux 
results in Sif being proportional to Y^^"^^^\ 

In summary, we have found that the scattering of two 
uncorrelated anyonic beams gives rise to HBT correla- 
tions which are bosonic in nature. This has been shown 
for a HBT interferometer threaded by an Aharonov- 
Bohm flux, and has also been suggested by the analysis 
of our toy model. The amplitude of the flux dependent 
cross-current correlations is non-analytic in the rates of 
the elementary two-anyon processes. Generalizing our 
model to finite temperatures [eV ~ hgT) allows QPs 
to tunnel from the inner edges to the outer edges, but 
otherwise no quantitative changes are expected. The ex- 
tension to finite range interaction will introduce higher 
harmonics at the fiux dependence [25]. More interest- 
ing is the inclusion of multi-channel edges (going beyond 
Laughlin fractions) , and eventually the generalizations to 
QPs satisfying non-Abelian statistics. 
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SUPPLEMENTARY MATERIAL 

Here we present a derivation of Pq.([5| and Eq.Q, the 
one-particle and two-particle rates, respectively. Let \ipi) 
and \ipf) be two many-body eigenstates of the system in 
absence of tunneling (the tunneling Hamiltonian is Ht). 
Very generally the transition rate between them due to 
the tunneling Hamiltonian can be written as 
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\{^P,\f\^f)\H{Ef-E,). 



(7) 



(8) 



Let us first consider the case of one-particle rate. For 
the sake of concreteness, we consider here tunneling 
through QPG A, all the other single particle rates be- 
ing similar. In this case \^f) is obtained by remov- 
ing a QP from edge 1 and transferring it to edge 3. 
Since we are interested in the total transition rate, we 
sum over all possible initial and final states. Notice 
that each edge is kept at a finite chemical potential /i^ 
[i = {1,2,3,4}) and that the initial states are weighted 
by w, = Z-i(V',|exp[-/3(i?o -EiM^^^)]IV'^), with Z = 
Trexp [— /3(i7o ^ Tli l^i-^i)]- To the lowest order in the 
tunneling amplitude the transition rate W^^-^f^ is given 
by 



W. 



(9) 

Here the operator H^'^ annihilates a quasiparticle on 
edge 1 and creates it on edge 3. Expressing the tun- 
neling operators in the interaction representation (with 
respect to Hq), Eq.([9]) can be rewritten as 

if 

(10) 



Notice that in Eq.(lO) we can extend the sum over final 



states to a sum over a complete set of states and obtain 



W, 



/GO 
-OO 



— oo 
oo 



sm 



{~ivt + Ic) 



The integration leads to Eq.([5]). In order to obtain 
Eq.(ll) we note that the two point correlation function 
for an edge kept at finite temperature and finite chemical 
potential is (this is equivalent to what presented in the 
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main text) 



^ilJ,v{t — x/v)/hj^l/3 



I3hv 



1/3 



The sum over the final states may be changed to a sum 
over a complete set of states; we can rewrite the expres- 



e'^ "^' ' {e- e ' " ' ■'""'"')^=^ion as 



I3hv 



(-i(vt - x) + Ic) 



(12) 



+00 



dtY: 



Let us now consider the total rate of transferring two 
quasiparticles from the external to the internal edges. 
Since there are no contributions to such a rate from sec- 
ond and third order terms in the tunneling amplitudes 
Fs, we need to consider the fourth order, we thus have 



x(^,|i7|(0) 



Wi 



1 



E,-Ho- iO+ 



T^H^pitm 



This may be rewritten as 



W. 



I = 



X {ipf\ Ht 



E,-Ho- iO+ 
Ht 5{Ef - E,) . (13) 



+ CX3 



dt 



dti 



dt-2 



{H^{ti)H'?{Q)Hf{t)H^\t + 12)) (18) 



Notice that in this case the many body eigenstate IV'/) 
is obtained from |^/)j) by transferring two quasiparti- 
cles. Being interested only in the lowest contribution to 
the current-current correlation modulated by the mag- 
netic flux, we study the contributions proportional to 
ir^rgrcrpl. Wehave 



The other three contributions proportional to 
T\TbV*cTv clearly read: 



II 



+00 



dt 



dt. 



dt^ 



H^{t,)H^iO)H^\t)H^^t + t2)) , (19) 



Ei - Ho - iO- 



III 



dt 



dti 



dty 



X {H'^{t,)H^{0)Hf(t)H^'{t + t2)) , (20) 



1 



-4tl 



xS{Ef-E,) + c.c.. (14) 

The above contribution corresponds to the rate 
{2, ABCV,^ tot (j)) of table III the corresponding ampli- 
tudes are represented in Fig.lsl Indeed the operator iJ^^ 
{H^^) annihilates a quasiparticle on edge 1 (on edge 4) 
and then creates it on edge 2 (3) respectively; similar 
statements apply to the operators H^^ and H^^ . Let 
us consider one of the four contributions proportional to 
r;^rBr5rx, obtained from Eq.([l4]), 



IV 



+00 



dt 



dti 



dt-i 



X {H'?{ti)H^{0)H^\t)H^\t + t2)). (21) 
We thus obtain 



(2,.AeC-D,*to,0)) 



JJ+2 



(I + II + III + IV) + c.c. (22) 



Again using Eq.(12) we can write Eq.(18) as 

T — r* r„r* r^p^27rw($„b+j *)/$o „-iei/V(L4 + Li)/'uft 



/3hv 



dt 



dti 



1= — ^wMAH^ r-ff^lV-f) 



1 



H^j.^\^,)SiEf~E,). (15) 



E,- Hn+ ^0^ 



Once again, moving to the interaction representation one 
can rewrite the previous expression as 



+ 00 



dt^w,{il),\H^{Q) 



1 



X sm 



X sm 



X sm 



X sm 



TT r. , il t2. 

\iv(—t H 

Phv^ ^ 2 2 ' 

— ii;(-i + - + 
pnv 2 2 



t^. 
2 ' 



iLi + Ic] 
iLs + Q 
iL2 + Ic] 



(23) 



X {^f\Hf{t) 



H^HtM) (16) 



E,-Ho+ ^0^ 
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where t has been shifted by (ti — t2)/2. Remarkably, 
changing variables in the terms II, III and IV, yields ex- 
actly the missing sectors in the ti and t2 integrals of con- 
tribution I. We can then combine the four contributions 
into a single expression 

I^n+III+IV ~ p-^PgP*p^g-2irij/(*„6+j3>)/<I>0g-iei/V(L4+Li)/-uft 

xsin|-^H-t-| + |)-iLi+g| ' (24) 

This integral can be evaluated explicitly in the limit 
iMiVB » 1 and (ii + L4 - L2 - L3) « hv/iveV) 
to obtain 

^|2,^BCi7,$,otO-)) ^ 2 ir^r^rcri,! cos [27ru{^r + ^ab + j *o)/$o] 
167rV3/y3^5/3 

^ 2l/3r(2/3)a4/3;i7/3^4/3 

= n ir^rsrcr^,! cos [27ri/($r + ^ab+j *o)/*o] , (25) 
where we write: 

r^rsr^ri, = ir^rsrcrj,! exp(-i27ri/$r/$o) • 

Note that in the main text we take $r to be zero for 
simplicity. 



